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Abstract 


Spatial-temporal pattern modeling of forest cover changes provides valuable information for better 
understanding the change process and determining the effective factors of areas under change. In this 
study, it was tried to identify the factors affecting the reduction and destruction of forests in western 
Iran using the capabilities of modern technologies including remote sensing, and then implementing 
in the form of a suitable computational model as a mathematical model based on the behavior of 
nature. In order to investigate the reduction of forest cover in Sardasht city of West Azarbaijan 
province, satellite images of MSS, ETM+ and OLI for the years of 1977, 2000 and 2018 were used. 
The above images were preprocessed, processed, and classified into two categories of forest and 
non-forest. The logistic regression method was used to study the relationship between forest cover 
reduction and physiographic and human factors. In order to obtain the lands suitability map, a 
logistic regression relationship was established between the forest cover reduction map from 1977 to 
2000 and 2000 to 2018 and also factors affecting it. Finally, a simple spatial model was proposed 
that was able to predict the spatial distribution of forest degradation using logistic regression. The 
results showed that over the 41 years, about 33721 hectares of forests of Sardasht city have been 
damaged. According to the results, it was determined that from topographic variables, parameters of 
distance from the road and distance from the village had the most impact on the forest degradation 
rate. 
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جغرافيا و مخاطرات محيطى. سال دهم شمارة چهلم. زستان ۱۶۰۰ صص ۹۹-۹۲ 
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پیش‌بینی daz‏ تخريب جنگل با استفاده از سيستم اطلاعات جغرافیایی و مدل ركرسيون لجستیک 
(مطالعه موردی: شهرستان سردشت) 
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حکیده 


مدل‌سازی الگوی مكانى - زمانی در تغییرات پوشش جنگل اطلاعات ارزشمندی را جهت درک بهتر فرآیند 
تغییر تعيين عوامل مؤثر مناطق در معرض تغيير فراهم می‌آورد. در اين پژوهش سعی و تلاش مابر آن 
بوده تا با استفاده از قابلیت‌ها و توانایی‌های فناوری‌های روز ازجمله سنجش‌ازدور» سیستم اطلاعات 
جغرافیایی عوامل موثر بر کاهش و تخریب جنگل‌های غرب ایران را شناسایی کرده و سپس در قالب 
الگوی محاسباتی مناسب به‌صورت یک مدل ریاضی بر گرفته از رفتار طبیعت پیاده‌سازی کنیم. در اين 
پژوهش به‌منظور بررسی کاهش پوشش جنگلی شهرستان سردشت در استان آذربایجان غربی, از طریق 
تصاویر سنجنده OLI; ETM* MSS‏ مربوط به سال‌های ATOT‏ ۱۳۷۹ و ۱۳۹۷ استفاده گردید. 
تصاویر ذکر شده مورد پیش‌پردازش و پردازش قرار گرفتند و به دو طبقه جنگل و غیرجنگل طبقه‌بندی 
شدند. برای بررسی ارتباط IG‏ کاهش پوشش جنگل با عوامل فیزیوگرافی و انسانی از روش آماری 
رگرسیون لجستیک استفاده شد. برای به دست آوردن نقشه شایستگی تبدیل اراضی, رابطه رگرسیون 
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لجستیک بين نقشه کاهش پوشش جنگلی سال‌های ۱۳۵۲ تا ۱۳۷۹ و ۱۳۷۹ تا ۱۳۹۷ و همچنین عوامل 
مؤثر بر آن» برقرار شد. درنهایت» یک مدل مکانی ساده که توانایی پیش‌بینی پراکنش مکانی تخریب جنگل 
را با استفاده از رگرسیون لجستیک als‏ ارائه شد. ee‏ نشان داد که در طول É!‏ سال حدود \ YYVY‏ 
هکتار از سطح جنگل‌های شهرستان سردشت کاسته شده است. با توجه به نتایج مشخص شد که از 
متغيّرهاى توپ و گرافی» پارامترهای فاصله از جاده و فاصله از روستا بر میزان تخریب جنگل بیشترین تأثیر را 
داشته‌اند. 


۱- مقدمه 


تخریب و نابودی زیستگاه یکی از مهم‌ترین دلایل انقراض بسیاری از گونه‌ها بوده است (سگن" و همکاران, 
57 بنابراین حفاظت از گونه‌ها مستلزم حفاظت از زیستگاه‌های آن‌هاست. اکوسیستم‌های جنگلی امروزه با 
تهدیدات بیشتری نسبت به ساير اكوسيستمها مواجه‌اند به‌طوریکه بررسی‌های FAO‏ نشان دهنده کاهش سه 
درصدی مساحت جنگل‌های جهان در حدود ۲۰ ساله اخير است (کینن "و همكاران. ۲۰۱۵). از طرف دیگر پاسخ 
اکوسیستم‌های جنگلی به تغييرات اقلیمی به گونه‌ای است که انجام اقدامات سازشی و حفاظتی را ضروری می‌سازد 
Ji)‏ "و همکاران» ۲۰۰۲). علاوه بر این» برخى از كونههاى جنگلی نظير انواع سنجاب گونه‌هایی به‌شدت وابسته و 
متکی به اين اکوسیستم‌ها به‌عنوان زیستگاه هستند و هرگونه تغيير و تخریبی در اين زیستگاه‌ها می‌تواند زمينه نابودی 
آن‌ها را فراهم آورد (کیرووسکی و همکاران؛ ۲۰۱۳). با توجه به موارد بالا می‌توان چنین نتیجه‌گیری کرد که امروزه 
حفاظت از جنگل‌ها به‌عنوان یکی از مهم‌ترین اکوسیستم‌های طبیعی که زیستگاه وسیعی را نیز برای GL‏ گوناگونی 
از حیات وحش فراهم می‌آورد بیش‌ازپیش احساس می‌شود. از ۱۶ میلیون هکتار مساحت کل جنگل‌های pal‏ 
تقريباً ۶۱ درصد متعلق به جنگل‌های ناحیه رویشی زاگرس است (جزیره‌ای و ابراهیمی‌رستاقی» ۱۳۸۲). اين جنگل‌ها 
وسيعترين رويشكاه بلوط ايران و قديمىترين جنگل بلوط در جهان با قدمت ۰ ساله هستند كه ١١‏ استان و 
حدود یک سوم جمعیت. بيش از نيمى از دام‌هاء حدود ۳۰ درصد از كل نزولات. ٤١‏ درصد از آب‌های جارى و ۳۲ 
درصد از آب‌های سطحى کشور را به خود اختصاص داده است (ثاقبطالبى و همکاران ۱۳۸۳). 
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برای GES‏ و ارزیابی اين تغییرات» استفاده از فنون و ابزارهای سنجش ازدور RS)‏ و سیستم اطلاعات 
جغرافیایی (GIS)‏ با برخورداری از امکانات تحلیلی می‌توانند نقش بسزایی در تولید اطلاعات مکانی داشته باشند 
(سومن "و همکاران, ۲۰۱۲؛ s hie i‏ همکاران, ۲۰۱۷)؛ بنابراین, در سال‌های اخیر به دلیل دسترسی آسان به تصاویر 
ماهواره‌ای و قابلیت‌های نرم‌افزار سیستم اطلاعات جغرافیایی» مدل‌سازی تغییرات پوشش ALS‏ و پیش‌بینی تغيبرات 
آن در ote!‏ رایج شده است و به اين منظور روش‌های متفاوتی چون تجزیه وتحلیل زنجیره مارکف. شبکه عصبی 
مصنوعی و روش‌های آماری مثل رگرسیون لجستیک گسترش يافته است (اسلام "و همکاران؛ ۲۰۱۸). مدل 
رگرسیون لجستیک ازجمله مدل‌هایی است که در مطالعات بسیاری برای مدل‌سازی رشد شهری و تخریب جنگل به 
کار گرفته شده است. برای نمونه» اشنایدر و پانتیوس" (۲۰۰۱) تغیبرات کاربری زمین در آبخیز ايبسوئيج در ايالت 
ماساجوست را در Sb‏ زمانی ۱۹۹۱-۱۹۷۱ مطالعه کرده‌اند و روش رگرسیون لجستیک را جهت مدل‌سازی تغییرات 
کاربری زمين به کار كرفتهاند. آن ها متغیرهای فیزیکی و اجتماعی برای فرایند مدل‌سازی به کار گرفته‌اند و شرایط 
آیندۂ اراضی را پیش‌بینی کردند. نتایج کاهش اراضی جنگلی و توسعةٌ مناطق مسکونی جدید در i ils‏ مطالعاتی را 
طی زمان نشان داد. جنك و ماسر"(۲۰۰۳) روش رگرسیون لجستیک را برای مدل‌سازی الگوهای رشد شهری در 
شهر ووهان در جين به کار گرفتند. نتایج اين مطالعه نشان داد نزدیکی به جاده‌ها و مناطق مسکونی موجود نقشی 
مهم در گسترش شهرها دارد. آن‌ها در اين مطالعه بیان کردند که رشد شهرها نقشی مهم در کاهش اراضی کشاورزی 
دارد. همجنين؛ لین“ و همکاران (۲۰۰۷) تغييرات کاربری زمین را با روش رگرسیون لجستیک در یک حوضة آبخیز 
در شمال تایوان مدل‌سازی کردند. مقادیر ROC‏ برای مدل‌سازی تغییرات هر یک از کاربری‌های کشاورزی» شهر 
اراضی جنگلی و مرتعی به ترتیب» ۰/۷۳۳ ۰۰/۹۸۳ ۰/۸۸ و ۰/۷۵۷ به دست آمد. اين مقادیر نشان داد مدل رگرسیون 
لجستیک توانایی مدل‌سازی تغییرات کاربری زمین را در منطقۀ مطالعاتی دارد. میلینگتون" و همکاران (۲۰۰۷) نیز در 
مطالعه‌ای فنون رگرسیونی را برای مدل‌سازی تغیبرات کاربری زمین در یک آبخیز مدیترانه‌ای به کار گرفتند. آن‌ها 
تصاویر ماهواره‌ای را جهت تهیه نقشه‌های کاربری زمین و روش مقایسه يس از طبقه‌بندی را برای آشکارسازی 
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رگرسیون لجستیک به خوبی قادر به تعيين تغییرات ایجاد شده در کاربری زمين است. از ple‏ مطالعات انجام شده در 
اين زمینه در حارج از کشور می‌توان به بررسی تغییرات پوشش جنگلی بين سال‌های ۲۰۰۵ تا ۲۰۱۲ در جنگل‌های 
کشور لائوس با استفاده از رگرسیون لجستیک توسط فوم‌فیلا و همکاران (۲۰۱۷) اشاره نمود. نتایج اين مطالعه نشان 
داد که تخریب جنگل رابطه معنی‌داری با ارتفاع از سطح درياء فاصله از جاده‌های اصلی و تغییر شیوه‌های کاشت 
دارد. در مطالعه‌ای دیگر سومن و همکاران (۲۰۱۰) به بررسی عوامل موثر بر عامل تخریب جنگل در جنگل‌های 
منطقه مرکزی میانمار پرداختند. در اين تحقیق نقشه تراکم تاج پوشش جنگل از سال ۱۹۸۹ تا ۲۰۰۷ با استفاده از 
تصاوير ماهواره‌ای به دست آمد و با استفاده از مدل رگرسیون لجستیک به مدل‌سازی مکانی تخریب جنگل با استفاده 
از fol ge‏ محیطی. پرداخته شد. نتایج OUS‏ داد ارتفاع و فاصله تا مناطق شهرىء بهره‌برداری و فاصله تا روستا به‌شدت 
بر احتمال تخریب جنگل می‌افزاید. از دیگر مطالعات انجام شده در این زمینه نیز می‌توان به مطالعات ویو" و 
همکاران (THE)‏ در کشور ویتنام و مطالعه ماهینی و ترثر "در حوضه آبخیز رودخانه بوروا (واقع در کشور استرالیا) 
اشاره نمود. آدهیکاری "و همکاران (۲۰۱۷) به بررسی تغیبرات کاربری اراضی و پوشش زمین منطقه کارناتکا در 
کشور هند پرداختند. بدين منظور از مدل رگرسیون لجستیک استفاده كردند. نتایج نشان داد که رابطه منفی بين 
تخریب جنگل و مناطق حفاظت شده جنگلی وجود دارد. 
در داخل کشور پیرباوقار *(۲۰۱۵) به مدل‌سازی تخریب جنگل با استفاده از رگرسیون لجستیک و GIS‏ در 
جنگل‌های غرب گیلان پرداختند. نتایج اين بررسی نشان داد که احتمال تخریب جنگل با متغیرهای شیب فاصله از 
جاده و مناطق مسکونی رابطه معکوس دارد. نتایج ارزیابی نشان داد که مدل به دست آمده با ضریب ROC‏ برابر 
۷ از تطبیق خوبی با واقعیت برخوردار است. استفاده ترکیبی از رگرسیون لجستیک و فرآیند تحلیل سلسله مراتببی 
(AHP)‏ نيز در مطالعه قجر* و همکاران (۲۰۱۲) به‌منظور برآورد سهم صخره و سنگ در حاکبرداری جاده‌های 
جنگلی در جنگل‌های آموزشی دانشگاه تربیت مدرس به انجام رسیده است. به اين منظون واحدهای زمین‌شناسی 
به‌وسیله فرایند تحلیل سلسله مراتبی (AHP)‏ ارزش‌دهی به همراه عامل شيب وارد مدل مربوطه شدند. از بين توابع 
اتصال مورداستفاده در قدرت برازش مدلء بهترین نتيجه مربوط به تابع اتصال Probit‏ بود. مهدوی و همکارانشان 
(۱۳۹۷) به مدل‌سازی تخریب جنگل‌های زاگرس با استفاده از رگرسیون لجستیک (مطالعه موردی: جنگل‌های 
چرداول استان ایلام) پرداختند. نتایج پژوهش آن‌ها OLS‏ داد که طی ۲۷ سال حدود ۱۰۳۳۲ هكتار از جنگل‌های 
شهرستان چرداول تخریب شده است؛ همچنین نتایح مدل‌سازی آن‌ها OUS‏ داد که متغیر جهت دامنه با دارا بودن 
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بيشترين ضريب تأثير CVT)‏ شايد مهمترين عامل بيوفيزيكى تأثي ركذار بر تخريب جنكل در منطقه موردمطالعه 
بوده است. پس‌ازآن, به ترتيب متغيرهاى شيب و ارتفاع از سطح دريا در تخريب احتمالى جنگل تأثير داشته است. 
شادمانی و همكاران (۱۳۹۹) به مدل‌سازی درجات تخريب جنگل‌های حوضه ماسال استان كيلان با استفاده از 
ركرسيون لجستيك يرداختند. نتايج مدل OLS‏ داد كه مدل ساخته شده با تابع اتصال Probit‏ از قابليت مناسبی در 
مدل‌سازی درجات تخریب برخوردار است. همچنین نتایج نشان داد که ۲۰۷۲۰ هكتار از جنگل‌های حوضه 
موردمطالعه با درجات بسیار کم تا متوسط با تخریب مواجه شده‌اند. 

در اين پژوهش سعی و تلاش ما بر آن بوده تا با استفاده از قابلیت‌ها و توانایی‌های فناوری‌های روز ازجمله 
سنجش ازدور» سیستم اطلاعات جغرافیایی عوامل موثر بر کاهش و تخریب جنگل‌های غرب ايران را شناسایی کرده 
و سپس در قالب الگوی محاسباتی مناسب به‌صورت یک مدل ریاضی بر گرفته از رفتار طبیعت پیاده‌سازی کنیم. 
درنهایت به شمای کلی و نزدیک به واقعیت از فرآیند تغییر در جنگل دست el‏ به‌طوری‌که بتوان از آن به‌عنوان 


راهنما برای انتخاب بهترین راه‌حل و ارائه راهکارهای مناسب جهت bie‏ و حراست از جنگل استفاده کرد. 


۲- مواد و روش‌ها 

۱-۲- داده‌های مورداستفاده 

در اين پژوهش از داده‌های سنجش‌ازدور شامل تصاویر MSS‏ (سال ۱۳۵۳/۶/۳۱ 'ETM*‏ (۱۳۷۹/۵/۱۱۲) و 
۱۳۹۷/۵/۱۷(۳۵]7) از سایت Earth explorer‏ و نيز نقشه‌های تویوگرافی ۱:۵۰۰۰۰ و عکس‌های هوایی ۱:۲۰۰۰۰ 
منطقه تهیه شده توسط سازمان جغرافیایی نیروهای مسلح استفاده شد. زمان انتخاب تصاویر» فصل اوج رشد گیاهان 
eal‏ 

-Y-Y‏ منطقه موردمطالعه 

منطقه موردمطالعه در جنگل‌های زاگرس شمالی و در استان آذربایجان‌غربی واقع شده است. منطقه موردمطالعه با 
موقعیت جغرافیایی ۳۷۳۹ عرض شمالی و LOOT‏ طول شرقی شهرستان سردشت است که جنگل‌های أن بين 1۰ تا 
۰ هزار هکتار تخمین زده می‌شود. شکل ۱ موقعیت منطقه را نشان می‌دهد. منطقه موردمطالعه در سیستم شبکه 
مسطح UTM’‏ در زون ۳۸ شمالی واقع شده است. 


1 Multispectral Scanner 

2 Enhanced Thematic Mapper 
3 Operational Land Imager 

4 Universal Transfer Mercator 


vé‏ جغرافیا و مخاطرات محیطی شمارة چهلم 
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شکل ۲- موقعیت شهرستان سردشت در کشور و استان آذربایجان‌غربی 


YY‏ روش تحقیق 

۱-۳-۲ تصحیحات هندسی و رادیومتریک 

جهت آماده‌سازی تصاوین ابتدا تصحیح‌های هندسی و رادیومتریک بر روی تصاوير انجام گرفت. در ابتدا با 
استفاده از نقاط کنترل زمینی )06 نقطه) از نقشه‌های توپوگرافی و استفاده از سیستم موقعیت‌یاب جهانی (GPS)‏ اقدام 
به تصحیح هندسی تصاویر گردید. خطاهای به‌دست‌آمده به ترتیب برای تصاویر اول» دوم و سوم برابر با ۸/۳۶ ۰/۱۷ 
و ۰/۵۱ پیکسل بود که بسیار مطلوب است. در مرحله بعد با کمک روش نمونه‌برداری مجدد" به روش نزدیک‌ترین 
همسايه ' تصاویر در جهت شمال توجیه و گردانده شدند. تصحیح رادیومتریک زمانی انجام می‌گیرد که از تصاویر 
چند “ale;‏ یعنی تصاویری که مربوط به فصول يا سال‌های مختلف و يا سنجنده‌های مختلف هستند. استفاده شود. 
دو نوع تصحیح رادیومتریک وجود دارده تصحیح رادیومتریک مطلق و تصحيح رادیومتریک نسبی. روش تصحیح 
رادیومتریک glee‏ نیازمند ورود داده‌های مربوط به خصوصیات اتمسفریک و کالیبراسیون سنجنده است. انجام اين 
تصحیح در اغلب موارد و بهخصوص برای داده‌های قدیمی کار بسیار سختی است (دندونکر * و همکاران» ۲۰۰۹) در 


مقابل تصحیح رادیومتریک نسبی با هدف کاهش متغیرهای اتمسفریک و غیره منتظره در ميان تصاوير چندزمانه 


1 Resampling 

2 Nearest neighborhood 
3 Multi Temporal images 
4 Dendoncker 
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صورت می گیرد. یکی از روش‌های تصحیح رادیومتریک نسبی کاهش تیرگی يديده' است. در حالت ایده‌آل پدیده- 
های تیره‌رنگ دارای تابش صفر در همه طول‌موج‌ها هستند. در اين روش فرض می‌شود كه در هر باند از تصوير 
می‌توان پیکسل‌هایی يافت که مقادیر آن‌ها صفر يا نزدیک به یک است (مثل (OT‏ به‌این‌ترتیب اثر اتمسفری تابش 
انحرافی به‌صورت یک مقدار ثابت به پیکسل‌ها در هر باند اضافه می‌گردد. به همین دلیل جهت حذف خطای 
رادیومتریک بایستی ارزش پیکسل‌های هر باند از حداقل by DN‏ به هر باند کم شود. کاهش تیرگی پدیده یک 
روش ساده است که به‌طور گسترده در بسیاری از موارد به کار گرفته می‌شود (چاوز " 2997). در اين مطالعه از اين 
روش جهت تصحیح رادیومتریک تصاویر موردنظر استفاده شده است که به‌راحتی در نرم‌افزار ENVI‏ قابلاجرا 
است. 

۲-۳-۲- انتخاب نمونه‌های تعلیمی 

برای طبقه‌بندی تصاویر ماهواره‌ای سال‌های مربوطه از روش طبقه‌بندی نظارت شده استفاده شد. در این پژوهش 
نمونه‌های تعلیمی به‌عنوان الگوی مشخصات طیفی طبقات. بر اساس کار میدانی انتخاب شدند. با توجه به هدف 
مطالعه در طبقه‌بندی (جنگل و غیرجنگل) نمونه‌های تعلیمی متناسب با پوشش هرکدام از کلاس‌ها در منطقه به 
تعداد لازم و پراکنش متناسب در سطح منطقه انتخاب شده‌اند. انتخاب نمونه‌های تعلیمی با توجه به پوشش جنگلی 
مناسب منطقه و قابل تشخیص بودن أن از مناطق غیرجنگلی با تفسیر بصری تصوير ترکیب رنگی كاذب و با استفاده 
از نقشه‌های تویوگرافی با مقیاس ۱:۵۰۰۰۰, نمونه‌های تعلیمی برای هر كلاس تعریف شدند. دو كلاس کاربری» 
به‌صورت جنگل و غیرجنگل در منطقه مشخص شد. سپس نمونه‌های تعلیمی به شیوه‌ی رقومی كردن روی صفحه 
تولید و تفکیک شدند. 

برای انجام طبقهبندى تصاویر از سنجنده اماساس و ای‌تی‌ام‌پلاس و OLI‏ مربوط به سال‌های ATOT‏ ۱۳۷۹ و 
۷ از روش‌های طبقه‌بندی نظارت‌شده استفاده گردید. يس از انجام دادن طبقه‌بندی تصاویر هر سال با روش‌های 
طبقه‌بندی نظارت‌شده تصوير با صحت بالاتر (روش ماشین بردار پشتیبان) جهت انجام عملیات‌های بعدی انتخاب و 
مورداستفاده قرار گرفت. در روش طبقه‌بندی نظارت شده محقق بر اساس واقعیت زمینی نمونه‌های تعلیمی را 
انتخاب و طبقه‌بندی را انجام می‌دهد. 


1 Dark Object Subtraction 
2 Chavez 
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۳-۳-۲- نهيه نقشه تغییرات و تخريب سطح جنگل مربوط به سالهاى مورد JB‏ 

پس از طبقهبندى تصاوير و ارزيابى صحت نتایج دقيقترين نقشههاى جنگل و غير جنكل مربوط به سالهاى 
WAV ۰۱۳۷۹ ۲7‏ منطقه موردمطالعه انتخاب شدند (بر اساس مقايسه صحت روشهاى مختلف طبقه‌بندی). 
نقشه‌های مذکور به‌صورت دو دوره از ۱۳7 تا ۱۳۷۹ و ۱۳۷۹ تا ۱۳۹۷ در محیط GIS‏ تقابل ۱ داده شدند. نقشه‌های 
pole‏ از عملیات تقایل نقشه تغیبرات پوشش جنگلی در دوره‌های ۱۳۷۹-۱۳۵۳ و ۱۳۹۷-۱۳۷۹ هستند که میزان و 
مکان تخریب پوشش جنگلی در دوره مورد نظر را OUS‏ می‌دهند. 

۳-۲--تهیه نقشه‌های فاکتورهای مؤثر در تخریب جنگل 

شناخت و جمع‌آوری اطلاعات مربوط به عواملی که در وقوع تخریب جنگل نقش مؤثرى دارند اولین مرحله 
مطالعه مدل‌سازی تخریب جنگل است. در اين پژوهش تعیین اين عوامل از طریق بررسی نتایج تحقیقات گذشته در 
مناطق مشابه صورت گرفت. يس از بررسی‌های انجام‌شده در منطقه موردمطالعه» اين عوامل در تخریب جنگل‌های 
منطقه مهم‌تر و موثرتر شناخته شدند که شامل: ارتفاع جهت شیب فاصله تا آبراهه شيب زمین, فاصله تا روستاء 
فاصله تا جاده می‌باشند. در اين راستا نقشه عوامل فوق با استفاده از نرم‌افزارهای جی‌آی‌اس ۱۰,۶ و انوی ۵۳ و 
ایدریسی سلوا تهیه شده و جهت تجزیه‌وتحلیل‌ها مورداستفاده قرار گرفتند. 

برای بررسی عوامل تأثیرگذار در تخریب جنگل‌های موردمطالعه. نقشه‌های مشخصه‌های ارتفاع از سطح درياء 
جهت شیب. فاصله تا آبراهه» شیب فاصله تا روستاها و جاده‌ها می‌باشند که در محیط ag GIS‏ شدند. با توجه به 
تعداد و پراکندگی خاص مناطق مسکونی و جاده‌ها و آبراهه‌های موجود در سطح منطقه تعداد ٠‏ حريم ۲۰۰ متری 
از مناطق مسکونی و ۱۰ حریم ۲۰۰ متری از جاده‌ها و ۱۰ حریم ۱۰۰ متری از آبراهه‌ها ایجاد شدند. 

0-۳-۲- تجزیه‌وتحلیل داده‌ها با رگرسیون لجستیک 

کاربرد رگرسیون لجستیک تعيين احتمال وقوع هر يك از سطوح متغیر کیفی دو حالته بر اساس متغيرهاى مستقل 
است و شرايط استفاده از آن به گونه‌ای است كه متغير وابسته حتماً يك متغيّر دو حالته است. ركرسيون لجستیک 
روش آماری متعلق به كروه مدلهاى آمارى خطى تعميم يافته است كه با استفاده از متغيّرهاى مستقل» احتمال وقوع 
یک حادثه را مورد پیش‌بینی قرار مىدهد. نكته مهم و اصلى در ركرسيون لجستيك اين است كه متغير وابسته یک 
متغير دوحالته است. يعنى اينكه مى تواند فقط عدد صفر به معنى عدم وقوع و يك به معنى وقوع رویداد باشد. در 
مورد تهيه نقشه حساسيت يا احتمال تخريب جنگل هدف ركرسيون لجستيك يافتن بهترين مدل برای تشريح روابط 


ميان حضور و يا عدم حضور متغير وابسته (تخريب جنگل) و مجموعه گروه‌هایی از متغير مستقل است. ركرسيون 


1 0 
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لجستیک از روش برآورد حداکثر احتمال ! (MLE)‏ برای بيدا كردن بهترین مجموعه پارامترهایی که مدل را بهتر 
برازش می‌کنند. استفاده می‌کند. خروجی مدل. ضریب‌های بين ٠‏ و ۱ خواهد داشت که از طریق تثورى فازی به 
احتمالات بالاتر از ۰/۵ ارزش Y‏ (تخریب) و کمتر از ۰/۵ ارزش صفر (بدون تخریب) می‌دهد و نقشه بولین تخریب 
را تولید می‌کند. رگرسیون با اين فرض بكار می‌رود که احتمال یک بودن متغيّر وابسته از منحنی لگاریتمی پیروی 
XS p‏ و مقدار آن توسط رابطه (۱) تخمین زده می‌شود (ایستمن F‏ ۲۰۱۲): 


exp(EBX) (\)‏ 1 20اه 


P(y = 1IX) = g 

AS پارامترهای برآورده شده و ۷ همان متغيّر وابسته است‎ B متغيّر مستقل»‎ X احتمال یک بودن متغیر وابسته»‎ P 
در اینجا تخریب جنگل مد نظر است.‎ 

به‌منظور حطی‌سازی رابطه Mb‏ تغییر لگاریتمی روی Ul‏ صورت می‌گیرد و به‌صورت Y daly‏ حاصل می‌شود. 

Loge (P/(1+P) = bo +b1x1 +b2x2 +... 00 + error term (Y) 

اين تغيير لگاریتمی سبب می‌شود احتمال پیش‌بینی شده در دامنه ٠‏ تا ۱ پیوسته باشد و خروجی مدل به‌صورت 
یک نقشه پیش‌بینی مکانی احتمال تخریب ارائه شود. 

مدل به دست آمده با استفاده از درصدی از نقاط بر گرفته شده از نقشه اولیه تخریب به روش‌های نمونه‌برداری 
سیستماتیک يا تصادفی طبقه‌بندی شده به‌صورت ارائه آماره‌های ROC‏ و Pseudo-R2‏ (مشخصه نسبی اجرایی) 
ارزیابی می‌شود. Pseudo-R2‏ برابر ۰/۲ به‌عنوان برازش مناسب در نظر گرفته می‌شود. ROC‏ آماره مناسبی برای 
ارزیابی اعتبار مدل است و می‌توان از OF‏ برای مقایسه تصوير مدل‌سازی شده با تصوير واقعی تخریب استفاده کرد. 
ارزش ۱ نشان دهنده توافق مکانی کامل و ارزش ۰/۵ نشان دهنده توافق کم مدل با واقعیت است (کلارک "و 
همکاران» (WI‏ 

از ابزارهایی که می‌تواند مورد استفاده برنامه‌ریزان در جهت کنترل روند تغییرات پوشش جنگلی قرار كيرد مدل‌ها 
هستند. برای برآورد پراکنش کاهش پوشش جنگلی شهرستان سردشت و تعين عوامل مؤثر بر آن از روش آماری 
رگرسیون لجستیک استفاده شد. رگرسیون لجستیک نسبت به ساير روش‌های آماری چند متغیره دارای مزایایی است. 
این مدل علاوه بر اينكه تابع» ارتباط مناسبی را با مدل رگرسیون خطى معمولی برقرار می‌کند متغیرها در اين روش 
می‌توانند هم شامل متغیرهای پیوسته و هم متغیرهای گسسته و يا هرگونه ترکیسی از اين دو نوع متغیر می‌باشند. 
همچنین اين مدل نسبت به دیگر روش‌های آماری چند متغیره به فرضیات کمتری نیاز دارد (لی و سامپته 4 ۲۰۰5). 


1 Maximum Likelihood Estimation 
2 Eastman 

3 Clark 

4 Lee 
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۳- نتایج و بحث 

برای آگاهی از وجود يا عدم وجود خطاهای هندسی در تصاویر از خطوط رقومی شده مربوط به آبراههها و 
جاده‌های موجود در منطقه کمک گرفته شد و خطای هندسی مشاهده نشد. همچنین به‌منظور کنترل دقت هندسی 
تصاویر حطوط رقومی مربوط به جاده‌ها از نقشه‌های رقومی مربوط به توپوگرافی ۱:۰۰۰۰ بر روی تصاویر تطابق 
يافته انداخته شد و تطبیق مسیر راه‌ها در روی تصوير با حطوط وکتوری۱ مورد مقایسه قرار كرفت و دقت تطابق 
خوب ارزیابی گردید. 

نتایج ارزيابى صحت نقشه‌های طبقه‌بندی در جدول ۵ ۲ و IY‏ شده است. نتایج طبقه‌بندی نظارت شده در 
شهرستان سردشت از نظر صحت طبقه‌بندی با استفاده از ضرایب صحت IS‏ و ضريب كاياء مورد مقایسه و آنالیز 
آماری قرار گرفت. به‌طوری که تصاویر سال‌های ۰۱۳۵7 ۱۳۷۹ و ۱۳۹۷ به ترتیب دارای صحت کلی AVM‏ ۹۹/۸۷ و 
۸ درصد است. با توجه به مقادیر ۹٩ CYL‏ درصد صحت كلىء می‌توان از اين نقشه‌ها برای تهیه نقشه تخریب 
استفاده کرد. برای ارزیابی صحت تصوير سال ۱۳۵۹ بخشی از منطقه که در هر سه دوره تغییر نکرده انتخاب شده و 
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شکل ۱۲-نقشه طبقات فاصله از آپراهه‌ها در شهرستان سردشت 
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به‌منظور ارزیابی مدل تهیه شده از شاحص ویژگی عملیات نسبی (ROC!)‏ و ضریب تشخیص کاذب " استفاده 


شده است كه به ترتيب براى دوره اول و دوم E‏ ۰۸۵ ليتق و ۱/۲۸۵۹ است. به دلیل قرار داشتن در 
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صحت مدل بيشتر مىشود و براى ضريب تشخيص كاذب عددى بين ۰/٤-۲‏ است كه تأییدکننده قابليت نسيتاً 


خوب مدل است. 


جدول 7- ارزيابى صحت مدل رگرسیون لجستيك با استفاده از شاخص‌های Pseudo-R? ; ROC‏ 
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فاصله تا آبراهه«۰/۰۰۰۰۱- شیب زمينءا ۰/۰۱۹۶- فاصله تا جاده× ۰/۰۰۰۶- فاصله تا مناطق مسکونی × 


وت 


ci‏ نتيجه گیری 

مدل‌سازی مكانى یک ابزار مناسب برای درک بهتر علتهاى تغييرات كاربرى اراضى و پوشش زمینی محسوب 
می‌شوند. اين تحقيق با هدف پیش‌بینی پراکنش مكانى تخريب جنگل‌های زاگرس و به‌منظور رديابى عوامل مؤثر بر 
تخريب جنكل صورت كرفت. 

در اين پژوهش نقشه تغييرات پوشش جنگلی مربوط به دو دوره ۱۳۵۳ تا ۱۳۷۹ و ۱۳۷۹تا WAV‏ به‌عنوان متغير 
وابسته و داده‌های رقومی ارتفاع شيب زمين» جهت شيب (ue)‏ فاصله تا coslo‏ فاصله تا روستاها و فاصله تا آبراهه 
به‌عنوان پارامترهای مؤثر در روند کاهش پوشش جنگلی به‌عنوان متغیرهای مستقل در برقراری رابطه رگرسیون 
لجستیک بكار گرفته شدند. خروجی رگرسیون لجستیک با شاخص Pseudo-R2‏ برای دوره ۱۳۵۹ تا ۱۳۷۹ برابر با 
۲ و شاخص ROC‏ برابر ۰/۸۱۲۶ و برای دوره ۱۳۷۹ تا ۱۳۹۷ با شاخص Pseudo-R2‏ برابر با ۰/۲۸۵۹ و 
شاخص ROC‏ برابر با ۰/۸۵۲۰ نشان‌دهنده برازش خوب مدل‌های به‌دست‌آمده با کاهش پوشش جنگلی واقعی و 
توانایی مناسب مدل در برآورد تغییرات جنگل در منطقه است. نتیجه حاصل شده با نتيجه برازش باقری و 
شتایی(۱۳۸۹) مطابقت دارد. 

با توجه به ضرایب منفی متغیرهای فاصله از مناطق مسکونی و جاده‌ها می‌توان كفت با کمتر شدن فاصله نسبت 
به اين عوامل توسعه مناطق انسان‌ساخت. جاده‌سازی و افزايش جمعیت نقش مؤثرترى در تخریب جنگل‌های منطقه 
موردمطالعه ایفا می‌کنند. ماتيو و همکاران QE)‏ گرونبرگ (۲۰۰۰) و ميريام و تیلور Y)‏ جاده‌سازی را به‌عنوان 
عامل اصلی تخریب جنگل در مناطق موردهطالعه خود معرفی کردند. یکی دیگر از عوامل مستقلی که در اين تحقيق 
موردبررسی قرار كرفت و در مدل به آن ضريب تعلق گرفت. فاکتور شيب است. در اين بررسی در مدل به فاکتور 
شیب ضریب منفی تعلق گرفت. که نشان دهنده ارتباط معنی‌دار ميان تغییرات شيب و تخریب جنگل است یعنی با 
افزايش شیب از میزان تخریب جنگل‌های منطقه کاسته می‌شود. به دلیل اينكه تخریب عمده جنگل‌های منطقه در اثر 
فعالیت‌های اگروفارستری اتفاق افتاده است و شیب‌های پائین برای کشاورزی مناسب می‌باشند و درنتیجه مناطق 
مذکور بیشتر مورد تخریب واقع شده‌اند. در ارتباط با تأثير عامل ارتفاع از سطح bo»‏ روی تخریب جنگل‌های منطقه 
موردمطالعه نتایج نشان داد كه در ارتفاعات پائین تخریب کمتر است که در خروجی رگرسیون لجستیک به آن 
ضریب مثبت تعلق گرفت؛ که نتيجه حاصله با نتایج تحقیقات پیرباوقار (۱۳۸۳) و رفیعیان (۱۳۸۵) که بیان کردند که 
در ارتفاعات بالا میزان تخریب جنگل بیشتر از ارتفاعات پائین است» هم خوانى دارد. در غرب به دلیل کوهستانی 


بودن در همه نقاط ارتفاعی بخصو Ue‏ یایین‌بند پراکنده است و تفاوت ناجيز اختلاف ارتفاعی بين مراكز جمعیتی و 
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جنگل وجود دارد. درحالی‌که در شمال کشور با افزايش ارتفاع به‌طور چشم‌گیری روستاها و مراكز جمعیتی کاهش 
می‌یابد. ضریب جهت شیب نیز مثبت بوده و طبق بررسی صورت گرفته از منطقه موردمطالعه. مشاهده شد که تراکم 
درختان جنگلی در جهت‌های شمالی بیشتر و وجود لکه‌های غیرجنگلی در این مناطق نیز به نسبت کمتر از ساير 
جهت‌ها است. جهت‌های شرقی به Wo‏ برخورداری از كرما و پوشش علفی بیشتر نسبت به جهت‌های شمالی در طول 
زمان بیشتر مورد تعرض قرار گرفته و جنگل‌های واقع در اين مناطق بیشتر تخریب یافته‌اند. 

با توجه به نتایج به‌دست آمده نرخ متوسط سالانه تخریب جنگل در شهرستان سردشت که مساحت کل آن 
۲ هکتار است. ۸۸۷/۳۹ هکتار در سال برای کل دوره زمانی ۶۱ ساله موردمطالعه برآورد می‌شود. با توجه 
به میزان برآورد شده می‌توان گفت» به‌طور متوسط سالانه ۰/۹۵ درصد سطح اولیه جنگل در منطقه کاسته شده است؛ 
بنابراین میزان تخریب صورت گرفته در اين منطقه بیشتر از متوسط جهانی (۰/۲ درصد) آن است. ازجمله دلایل اين 
کاهش ۳۳۷۲۱ هکتاری از طبقه جنگل در طی ۶۱ سال را می‌توان استفاده‌های نادرست. جاده‌سازی‌های زياد به دلیبل 
مناطق مرزی و جنگ‌های متعدد. استفاده نادرست تفرجی, تحولات جمعیت» گسترش شهری و کارهای عمرانی و 


اثرات ناشی از OF‏ در شهرستان سردشت اشاره کرد. 
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